Singularities of Hele-Shaw Flow and Shock Waves in Dispersive Media 
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We show that singularities developed in the Hele-Shaw problem have a structure identical to 
shock waves in dissipativeless dispersive media. We propose an experimental set-up where the cell is 
permeable to a non-viscous fluid and study continuation of the flow through singularities. We show 
that a singular flow in this, non-traditional cell is described by the Whitham equations identical to 
Gurevich-Pitaevski solution for a regularization of shock waves in Korteveg-de-Vriez equation. This 
solution describes regularization of singularities through creation of disconnected bubbles. 
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1. Introduction. A broad class of non-equilibrium 
growth processes in two dimensions are characterized by 
a common law: the velocity of the growing interface is 
determined by the gradient of a harmonic field (often re- 
ferred to as Laplacian growth) p|. 

Growth processes determined by a harmonic field, 
where no cut-off scale is introduced, are important. The 
theory of this kind of growth is deeply related to fun- 
damental aspects of conformal maps, integrable systems, 
2D quantum gravity and random matrices 2J. However, 
singularities makes this problem ill-defined and thus not 
achieved experimentally. 

The goal of this Letter is two-fold. One is to suggest an 
experimentally achievable set-up, where a high-rate flow 
can continue through singularities, without being curbed 
by dissipative forces. Another is to stress a deep and 
important parallel between the singularities of growing 
interface and "gradient catastrophes" known in the the- 
ory of shock waves in dispersive media. The Korteveg-de- 
Vriez (KdV) and other nonlinear waves equations feature 
this phenomenon. This relation allows one to effectively 
study complicated singular processes. 

The relation between the Hele-Shaw problem and dis- 
persive nonlinear waves suggests a unique regularization, 
namely one which plays a similar role as dispersion of 
nonlinear waves. There, dispersion, (no matter how small 
it is) becomes crucial near a singularity, converting shock 
waves to oscillatory structures 

2. The Hele-Shaw flow. The Hele-Shaw cell is a nar- 
row space between two plates filled by incompressible 
viscous liquid (say oil). Air (regarded as inviscid and 
incompressible) occupies a part of the cell forming one 
or several bubbles. Traditionally air is pumped into one 
bubble, while oil is extracted from the cell at a constant 
rate Q > through the edges placed at infinity. Without 
surface tension, the interface may develop singular cusps 
at a finite time pj . At this moment the problem becomes 
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FIG. 1: Schematic Hele-Shaw experimental set up. Shaded 
region represents a small bubble of air in the ambient oil. The 
insulator is permeable to air but not to oil. 



ill-posed. 

We suggest a novel set up where the upper plate is 
permeable to air and is connected to a reservoir of air. 
Neither plates are permeable to oil. This set-up may 
be achieved by sticking a Gortex-like material to glass 
plates having small perforations (FigGJ. We thus assume 
that all air bubbles are kept at the same atmospheric 
pressure. We also consider a retraction problem, where 
oil is extracted from the cell, i.e. Q < 0. 

In both cells the flow obeys the same local equations 
(D'Arcy's law): The local velocity of oil, averaged over 
the gap between the plates, is proportional to a gradi- 
ent of pressure p. In an incompressible fluid pressure is 
a harmonic function. Thus the normal velocity of the 
interface, v n , is proportional to the normal derivative of 
pressure, and in proper units v n — — d n p. If in addi- 
tion surface tension is ignored, the pressure is continuous 
across the interface, and may be set to p — within each 
bubble. Thus the pressure solves the exterior Dirichlct 
boundary value problem: Ap = 0, with p = on the 
boundary, and p — ► —Q log \z\ at infinity. 

3. Bubble break-off and merging. Let us now describe 
an air retraction process in our experimental set up. A 
typical local evolution of an air bubble consists of few 
phases illustrated in Fig. (i) As oil is injected, the air 
bubble contracts; (ii) The air bubble forms a singular nar- 
row neck and breaks-up into one or several disconnected 
bubbles; (iii) All bubbles subsequentially contract, while 
loosing air through the upper plate. 

As we show in this Letter, the method of regularization 
of shock waves provides a detailed description of retrac- 
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tion evolution in the limit of zero surface tension. This 
is in contrast with the injection process in a traditional 
Hele-Shaw cell where solutions at the zero surface ten- 
sion limit are valid only until a cusp is formed. This 
difference, suggests that surface tension is a non-singular 
perturbation in our new experimental setup. 

4- Dispersive regularization of shock waves. The term 
shock wave is usually used to describe nonlinear waves in 
presence of dissipation. In the absence of dissipation the 
shock- wave behavior is different. It is resolved by con- 
verting to oscillations. The method to study this compli- 
cated behavior was suggested by Gurevich and Pitaevskii 
(GP) 3], and refined in later works HHQ. They 
studied solutions to the KdV equation 

3 o 
®t 3 u — -ud tl u + ed tl u, (1) 

with step-like boundary conditions: u(ti — > — oo) = uq 
and u(t± — > +oo) = 0. Initially the function u{t^,ti) is 
smooth (ta = 1 of FigJSJj so the dispersion term df ' u 
may be neglected. However, the Hopf-Burgers equation 
d t3 u = |it<9 tl M, thus obtained, always develops a shock 
wave, a singularity with an infinite slope (£3 = of Fig. 
|3J) followed by an unphysical overhang (£3 = —1 of Fig. 
|2J). This signals that the limit e — > is singular. 

In fact the solution with finite e never develops a shock 
wave. Before the singularity occurs the wave breaks into 
fast oscillations (see Fig. OJ), of a period scaling with e. 

When e -C 1, the oscillatory regime can be described 
by slowly modulated periodic solutions of the KdV equa- 
tion. These are given by the elliptic function 

U (i 3 ,ii)~2a-dn 2 + Vt 3 ), m) +7, (2) 

which moduli a, m, V, 7 (i.e. frequency, amplitude etc.) 
depend on the times £3, t\. This dependence is described 
by the Whitham equations These equations appear 
below in (|11I12|> to describe the Hele-Shaw flow. 

A similar situation takes place in our problem. We 
will show that a typical flow is identical to averaging 
the solution over fast oscillations, and is, essentially, 
described by the same Whitham equations. 

5. Correspondence between interface dynamics and 
shock-wave solutions. In order to establish a link between 
modulated periodic solutions of the KdV equation and 
growth of planar domains let us recall the notion of the 
spectral curve fl5| . The spectral curve (or Riemann sur- 
face) encodes periodic solutions of nonlinear integrable 
equations. In the KdV case it is a hyperelliptic curve 

y 2 = R m (z), m = 2£+l, (3) 

where y and z are complex variables, and R m is a poly- 
nomial of an odd degree, such that its roots are real. The 
spectral curve of solution @ is given by a polynomial of 
degree m = 5 and eq. (|10|l . A non-periodic solution at 
small e is described by a slowly varying spectral curves. 



A shock wave indicates a singular evolution when the 
curve changes its genus. We will show that for the inter- 
face dynamics an increase of the genus implies a bubble 
break-off, since the interface is a real section of the curve 
when the coordinates z and y are real|8|-|l0j. 

6. The critical flow. Once injecting air into Hele-Show 
cell, the bubble develops a "finger". Its tip is pushed 
away with increasing velocity that eventually may result 
in a cusp-like singularity (Fig. [2J ■ 

Let us choose the origin at the point where a cusp 
would form, and simplify the argument assuming that 
the finger is symmetric with respect to its x-axis. In 
dimensionless units, where the size of the entire droplet 
is of order 1 we have \y\ <§; \x\ <C 1. Let us denote the 
distance between the tip and the origin by u(t) <C 1. It 
sets the only (time dependent) scale of the critical flow. 

A critical flow of an isolated finger is conveniently for- 
mulated as the time evolution of the Cauchy transform 
of air bubbles: 

h {±) (z) = ^- <f z = x + iy. (4) 

2m y 7 z — z' 

This integral defines analytic functions h (z) for z inside 
and outside air respectively. A direct calculation shows 
that D'Arcy law is equivalent to 

ftftW (z) = 0, fitfe(-) (z) = -2id zP (z, z) (5) 

The first equality implies an infinite series of conserved 
quantities first observed in [Tl|. 

A generic form of the critical finger (insets 7 and a 
in Fig. 0) is given by a curve ||, i/ 2 = R m (x), such 
that all the coefficients of the polynomial R m (x) are real 
and it has at least one real root x = u at the tip of the 
finger. In the critical regime the coefficients of the poly- 
nomial u m R m (x/u) are of order 1, while u — > as time 
approaches the critical time t c . There a finger becomes 
the (2, 21 + 1) cusp y 2 ~ \x\ 2e+1 . 

The critical flow is conveniently described in terms 
of the height function 0. The height function y(z) = 
\J R rn (z) is defined on a hyperelliptic Riemann surface, 
and is an analytic function outside the finger, having 
boundary value y(x) on the boundary of the finger. In 
it was proved that: (i) the finger remains self-similar, 
i.e., remains of the form of polynomial of a fixed degree, 
only if all its simple roots are real. In this case branch 
points other than u correspond to additional finite size 
droplets, if any. As a conseque nce, a sole finger is de- 
scribed by a degenerate curve \J R m (z) = v ' z ~ uPt(z), 
where Pi is a polynomial of degree t; (ii) The coefficients 
£3, . . . , £2^+3 (called deformation parameters) in front of 
all positive fractional powers of z in the expansion of 

«+i 1 

y{z) = y/RM) = 5) h k+1 z k -^+0(z- 3 / 2 ) (6) 

do not depend on time. Furthermore, the coefficient in 
front of the first negative power, z -1 / 2 , is proportional 



FIG. 2: Graph of numerical solution for branch points iti(ti) 
compared to the hodograph for u(ti). Finger shapes at early 
(7), singular (/3) and late (a) stages of evolution. 

to time t\ ~ Q(t — t c ). The negative tail of the series © 
depends on time in a non-trivial way. 

For simplicity let us consider a single finger and study 
the behavior of h ± (z) in the domain \u\ <C \z\ <C 1, far 
from the tip but around the finger, where details of the 
tip are not seen. The Cauchy integral J2J for z > u 
on the positive real axis contains only regular terms (i.e. 
positive integer powers of z). Therefore, M + ) extends 
analytically as a regular function to the whole domain of 
interest. The singular part of the function (contain- 
ing also fractional powers of z) then provides all necessary 
information. It has a cut inside the finger drawn along 
the :r-axis (see FigEJ. 

In order to compute the Cauchy integral (0J, we ex- 
pand the function y(x) in a series in half- integer powers 
of x and evaluate hS'[z) term by term. The function 
h(~\z) has a cut along the finger axis from u to infinity 
and takes opposite real values on the two sides of the 
cut. Conventionally, we choose a single- valued branch 
such that y(x + iO) > 0. For real negative z = x we find 
hf-~\x) — y{x). Therefore, the singular part of hS~\z) 
is to be identified with the height function 

h^\z)=y{z) = ^RjJ) (7) 

Now consider the evolution of the finger (or the curve 
©) encoded by Eq.JSJ. The leading behavior of pressure 
around the finger follows from solution to the Dirichlct 
boundary value problem around a slit (since \y\ -C \x\, 
the finger roughly looks like a cut): p(z, z) oc Re \J — z. 
This proves that higher terms in the expansion of the 
height function © are conserved. 

7. Dispersionless KdV hierarchy. We have reformu- 
lated the Hele-Shaw flow as an evolution of a hyperelliptic 
curve. Notably, the spectral curve of KdV equation 
evolves exactly in the same way. In particular, it follows 
from I15I7II that before the break-off the scale u evolves 
with time and other deformation parameters according 
to the dispersionless KdV hierarchy [t| 

2"«!9t 2 „ +1 w + (2n + \)\\u n d tl u = 0. 

8. (2,5) critical finger. A generic flow is represented by 
a family of critical fingers of the form J3J with 21 + 1 = 5 



FIG. 3: Solutions of the hodograph for tz = — 1 (green), 
(red) and 1 (black). Unphysical multivalued region (in dashed 
line) of the t$ = — 1-graph is replaced by an oscillatory solu- 
tion in heavy gray line (numerical). 



y 2 = (x-u)(x-d+) 2 (x-d-) 2 . (8) 

The a: -term, and therefore <5 in can always be elim- 
inated by a shift (2d + + 2d- + u = 0), leaving i 3 as the 
only deformation parameter. Substituting this to © we 
obtain the hodograph solution [l6l Il7| 

-u 3 + -ut 3 +h =0, (9) 

implicitly giving the branch point u in terms of t\, t%. 

A singularity occurs when u merges with one or two 
double points d± — j(— u ± \f— 24*3 — 5u 2 ). 

The features of finger's evolution crucially depend on 
the sign of £3 (Fig. 3). If £3 > 0, then the double points 
never reach the real axis. The function u(t\) is single- 
valued and in both extraction and injection processes the 
finger never becomes singular (inset a in Fig|2J). 

In the case t$ = the solution is u(t) oc — (Q(t— tc)) 1 / 3 . 
This corresponds to a finger which evolves into the (2, 5)- 
cusp y 2 = x 5 . At t = t c , all the three roots coincide: 
u = d+ = d- = 0. An interesting feature of the (2,5) 
cusp (and of all cusps (2, 2£ + 1) at even £) first noted 
in is that the evolution can be extended beyond the 
cusp by means of the same hodograph equation © ■ 

The most interesting case is £3 < 0. A formal solution 
in the injection case leads to the (2,3)-cusp, y 2 oc x 3 , 
which can not be continued. The plot in Fig. 3 becomes 
multi- valued in the region t\ < |(— i 3 ) 3 . 

9. Air extraction - bubble break-off. We treat air ex- 
traction, where Q < 0, and the physical time is —t%. At 
an early stage, the finger tip is far to the left (u is large 
negative and t\ is large positive). Extracting air from 
the bubble results in a tip motion from left to right and 
a changing shape of the finger in accordance with the 
hodograph solution. The evolution follows the branch 
of the plot from t\ = —00 until the point A in Fig. [2] 
At this stage, the double points d± of the curve © are 
complex. At t\ = tA = ^/54/5 (—t 3 ) 3 ^ 2 they approach 
the real axis and merge. As t\ — > tA, the finger devel- 
ops a thin neck around the point — it/4, which breaks at 
t\ = tA through the (2,4) cusp y 2 oc x A . At the break- 
off the curve is y 2 = (x — u)(x + ju) 4 (u is negative). 
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The bubble which breaks off from the main one has area 
2 u) 7 / 2 . After this singular point, the evolution ac- 
cording to the hodograph solution J^J is unphysical. It 
leads to an interface with a self-intersection point. 

The actual evolution at t\ < <a does not follow the 
cubic parabola Q in Fig. 3. The correct extension of the 
solution beyond the singularity describes a small bub- 
ble breaking off from the finger tip. A physical solution 
describes the interface (or a spectral curve) y 2 = -Rs(cc) 
with Rs(x) having only one rather than two double roots: 



This condition closes the system of equations (|f If) after 
the break-off. However, unlike the algebraic equations 
for tk, this condition is transcendental. A solution is 
available through elliptic functions. It gives the time de- 
pendence of functions u\, u 2 , u% plotted in Fig|2] After 
the break-off, the small bubble starts to evapo rate and 
eventually disappears at t\ = tc = — -\/2/27 (— t^] 3 / 2 
corresponding to the point C on the plot. After that the 
solution switches back to the cubic parabola. The finger 
proceeds further without obstacles. 



y 2 = (x - u\)(x - u 2 ){x - u 3 )(x ~ d ) 2 , (10) 

where do — V — ^(ui+u 2 + u 3 ). Now the height function 
y(z) has two cuts: one small cut inside the small bubble 
and an infinite cut inside the finger (inset (3 in Fig. 2). 
The double point do moves between them. Mathemat- 
ically, this means that the complex curve is of genus 1 
(an elliptic curve). Similarly to the genus-0 case JBJ, a 
substitution (|10fl into © we obtain 

12i x = Ul - 4£/3 , -12t 3 = Ul + 2U 2 , (11) 

where Uk = u\ + u\ + u 3 . 

In case of two bubbles, t\, t 3 are not enough to fix 
the dynamics. An additional conserved quantity fol- 
lows from the condition on pressure. Since air is under 
the same pressure in both bubbles, we write dp — 

2Re d z pdz = 0, or, using ©, Im dth^dz = 0. 
Since h^(ui) — y(ui) — 0, we rewrite this condition as 
dt ^Im yj R^{x) dx^j = 0. The quantity under the in- 
tegral is purely imaginary and vanishes at t\ = t^, hence 

f 3 VMx) dx = 0. (12) 



10. Discussion. Below we highlight our major results. 
We propose a novel, experimentally accessible modifi- 
cation of the Hele-Shaw cell, where the flow proceeds 
through singularities without being curbed by a surface 
tension. In this cell singularities of the flow are resolved 
by creation of new bubbles, all kept at the same constant 
pressure. This behavior can be seen experimentally. 



The flow in such an air-permeable cell can be studied 
with the help of mathematical tools of soliton theory. We 
identified bubble break-off with the shock wave behavior 
of the KdV equation. 
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